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ABSTRACT 

We report on the discovery of methyl acetate, CH3COOCH3, through the detection of a large number 
of rotational lines from each one of the spin states of the molecule: AA species (Ai or A2), EA species 
(Ei), AE species (E 2 ), EE species (E 3 or E4). We also report the detection, for the first time in space, 
of the gauche conformer of ethyl formate, CH3CH2OCOH, in the same source. The trans conformer 
is also detected for the first time outside the galactic center source SgrB2. From the derived velocity 
of the emission of methyl acetate we conclude that it arises mainly from the compact ridge region 
with a total column density of (4.2±0.5) x 10 15 cm -2 . The derived rotational temperature is 150 K. 
The column density for each conformer of ethyl formate, trans and gauche, is (4.5±1.0)xl0 14 cm~ 2 . 
Their abundance ratio indicates a kinetic temperature of 135 K for the emitting gas and suggests that 
gas phase reactions could participate efficiently in the formation of both conformers in addition to 
cold ice mantle reactions on the surface of dust grains. 

Subject headings: Molecular data — Line: identification — ISM: molecules — ISM: abundances — 
ISM: individual objects (Orion KL) 



1. INTRODUCTION 

The line survey of Orion carried out by Tercero and 
collaborators using the IRAM 30-m telescope presents 
a forest of lines arising from isotopologues and vibra- 
tionally excited states of abundant specie s (see, e.g., 
iTercero et aLllMLCl l20lH IDalv et al. ||20H . The prob- 
lem of identifying these features was a real challenge as 
initially we had more than 8000 unidentified lines. The 
analysis of the data has been done molecule by molecule 
(|Terceroll20T2] ). For each species we explored the lit- 
erature for spectroscopic information on the isotopo- 
logues and vibrationally excited states but substantial 
laboratory work was missing for most of these species. 
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In 2006 we started a close collaboration with differ- 
ent spectroscopic laboratories that allowed us to iden- 
tify nearly 4000 of these unknown lines (often called 
weeds). The 13 C and 15 N isotopolog ues from ethyl 
cyanid e, C H3CH2CN, were meas ured by iDemvk et al. I 
(120071) and lMargules et al. I(|2009f ). Several vibrationally 
excited states of t he same spec i es wer e characterized in 
the laboratory by IDalv et al. I (pOllh . The 13 C, 18 0, 
and deuterated isotopologues of methyl forma t e wer e 
observed in the labor atory by iCarvaial et alH (2009ft . 
iMargules et~aT~l (|2010f) . and ITercero et al. I (|2012f) ~The 
vibrationally excited state ^12 of formamide was mea- 
sured in the laboratory bv lMotivenko et al. I (2012). Fi- 
nally, Lopez et al. (in preparation) have characterized 
in the laboratory several vibrationally excited states of 
vinyl cyanide, CH2CHCN, which have been identified 
in Orion together with all its isotopologues 13 C, 15 N, 
and D. All these isotopologues and vibrationally excited 
states were detected in space for the first time thanks 
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to these new laboratory data. In addition, the study of 
the survey was divided in the analysis o f different fam- 
ilies of molecules: CS bearing s pecies (ITercero et al. I 
12010ft . silicon bearing molecules (ITercero et al. I l2011h . 
SO and S0 2 (|Esplugues et al. 1 12013 ) . Work on other 
species such as HCN, HNC, and HCO+, CH 3 CN, HC 3 N 
and HC5N, and organic saturated O rich species, is in 
progress (Marcelino et al., Bell et al., Esplugues et al., 
Lopez et al. in preparation). 

Although many of the still 4000 remaining U lines 
could belong to rare isotopologues of complex organic 
molecules, we think we are ready now to start the search 
for new molecular species. The study of a cloud such as 
Orion could provide important clues on the formation of 
complex organic molecules on the grain surfaces and/or 
in the gas phase. A systematic line survey with most 
weeds removed permits to address the problem of the 
abundances of isomers and derivates of key species, such 
as methyl formate, ethyl cyanide, and others. Moreover, 
it will constitute the best spectral template for future 
ALMA observations of hot cores. 

In this letter we report the detection for the first time 
in space of methyl acetate, CH3COOCH3, from the de- 
tection of many lines from the states AA, AE, EA, and 
EE (E3, and E4) of this molecule at 3, 2, and 1.3 mm. 
The gauche conformer of ethyl formate (an isomer of 
methyl acetate), for which t he anti conformer was pre- 
viously detected in SgrB2 bv lBelloche et al. I ()2009f ). has 
also been detected in Orion. This detection of ethyl 
formate outside the galactic center indicates that this 
species is also efficiently produced in hot cores. 

2. OBSERVATIONS 

The observations were carried out using the IRAM 
30-m radiote l escop e during several periods (see 
ITercero et al. I I2010D . System temperatures were in 
the range 200-800 K for the 1.3 mm receivers, 200-500 
K for the 2 mm receivers, and 100-350 K for the 3 
mm receivers, depending on the particular frequency, 
weather conditions, and source elevation. The intensity 
scale was calibrated using two absorbers at different 
temperatures and using the Atmospheric Transmission 
Model (ATM, Cernicharo 1985; Pardo et al. 2001). 
Pointing and focus were regularly checked on the nearby 
quasars 0420-014 and 0528+134. Observations were 
made in the balanced wobbler-switching mode. We 
pointed the observations towards IRc2 at a(J2000)=5' i 
35 m 14.5 s , <5(J2000)=-5° 22' 30.0". The data were pro- 
cessed using the IRAM GILDAS softwarfQ (developed 
by the Institute de Radioastronomie Millimetrique). 
We considered lines with intensities >0.02 K, covering 
three or more channels. Figures are shown in units 
of main beam antenna temperature, T ' mb = ^*a/vmb, 
where t\mb is the main beam efficiency. A more detailed 
descri ption of the observations can be found in iTercerol 
(I20L1 . 

3. FREQUENCY AND INTENSITY PREDICTIONS FOR 
METHYL ACETATE 

Met hyl acetate [CH3 - 0-C( =Q)- CH 3 ], has been stud- 
ied bv lSheridan et all (|1980l ) using Stark spectrometers 
in the frequency region from 8 to 40 GHz. In 2011, 

1 http: / /www.iram.fr /IRAMFR/GILDAS 
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Figure 1. Selected lines of Methyl Acetate at 3 mm, 
CH3COOCH3, towards Orion-IRc2. The lines from the different 
states are identified. The continuous green line corresponds to all 
lines already modeled in our previous papers (see text). 

the investigation was considerably improved with 315 
new lines recorded using the molecular beam Fourier 
transform microwave spectrometer in Aachen and 519 
lines recorded using the Free jet absorption Stark- 
modulatc d millimeterwave (FJA SMmm) spectrometer in 
Bologna (jTudorie et al. 1 120111 ) . A newly written pro- 
gram BE LGI-Cs-2Tops based on the Hamiltonian de- 
scribed bv lOhashi et al. I (|2004D was used to fit the com- 
plete data set with 27 molecular parameters, up to J 
= 19 and K a = 7. More than 800 new microwave 
and millimeter- wave measurements were assigned to the 
ground-state transitions in methyl acetate and fit, lead- 
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Figure 2. Selected lines of Methyl Acetate at 2 mm, 
CH3COOCH3, towards Orion- IRc2. The lines from the different 
states are identified. The continuous green line corresponds to all 
lines already modeled in our previous papers (see text). 

ing to root-mean-square deviations of 4 kHz for the mi- 
crowave lines and of 40 kHz for the millimeter-wave lines, 
i.e., to residuals essentially equal to the experimental 
measurement errors. The heights for two internal rota- 
tion barriers were determined to be of 102 cm -1 for the 
acetyl CH3 internal rotor and of 422 cm" 1 for the ester 
CH3. As the theoretical background and code presently 
used h ave been extensively described by iTudorie et al.l 
(|2011f) , we will summarize here only the main character- 
istics. 

The BELGI-Cs-2Tops program is restricted to i) asym- 
metric top molecules containing two non-equivalent CH3 
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Figure 3. Selected lines of Methyl Acetate at 1.3 mm, 
CH3COOCH3, towards Orion-IRc2. The lines from the different 
states are identified. The continuous green line corresponds to all 
lines already modeled in our previous papers (see text). 



internal rotors, and ii) molecules belonging to the Cs 
point group at equilibrium. This code is most closely 
related to a code used only once in the literature for a 
treatment of the microwave spectrum of t he molecule N- 
meth ylacetamide [CH 3 -NH-C(=0)-CH 3 ] (jOhashi et al. I 
2001 although some differences exist as detailed in 
Tudorie et all (poTTI) . 

The torsional and rotational Hamiltonian is diagonal- 
ized in separate symmetry blocks, each characterized by 
one of the five (<ri, 02) pairs, where o\ and a% desig- 
nates the symmetry indicators for each of the A and 
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Figure 4. Selected lines of the trans (red line) and gauche (blue) 
conformers of ethyl formate. CH3CH2OCOH, towards Orion-IRc2. 
The synthetic spectrum corresponds to the same column density for 
both conformers (see text). The continuus green line corresponds 
to all lines already modeled in our previous papers (see text). 

B tops. We follow the notation used in Table 1 of 
IQhashi et al~l (l200l . i.e. we have for the 5 symmetry 
species: ai—0, er 2 =0, AA species (Ai or A 2 in the Gig 
permutation-inversion group of the molecule), ai=±l, 
(72=0, EA species (Ei), oi=0, ct 2 =±1, AE species (E 2 ), 
<7i=±l, er 2 =Tl, EE species (E 3 ), cri=±l, ct 2 =±1, EE 
species (E4). The higher barrier hindering the ester 
methyl group internal rotation (422 cm -1 ) corresponds 
to the smaller splittings between the AA and the E 2 
lines, whereas the lower barrier (102 cm -1 ) hindering the 
acetyl methyl group is responsible for the larger splittings 
between the AA and the Ei lines. The splittings between 
the E 3 lines and the E4 lines are due to the coupling 



between the two tops. Statistical weights for methyl ac- 
etate are 16, 16, 16, 8 and 8 for the Ai o r A 2 (AA), Ei 
(AE) , E 2 (EA), E 3 , and E 4 , respectively (|Ohashi et al. I 
2004). The zero point energy for the J=K=0 levels are 
99.9450 cm" 1 , 99.9551 cm" 1 , 101.0947 cm" 1 , 101.1047 
cm -1 , and 101.1049 cm" 1 for A, E 2 , Ei, E 3 , and E 4 
species, respectively. 

The spectroscopic c onstants determined f rom the pre- 
vious fit published in iTudorie et al. I (|2011[ ) are used in 
the present study to predict the transition frequencies for 
the u t =0 torsional ground state up to J = 30. Since we 
did only fit the observed lines in the laboratory up to J — 
19, K a <7, large uncertainties exist for the higher J and 
K a values (we estimate la uncertainties of 0.1-0.5 MHz 
for J <25 and 0.5-2 for 25<J<30). For the present paper 
we also have implemented in the BELGI-Cs-2Tops the 
calculation of the intensities usin g the same method as 
described in lHougen et al. I (|1994t ): lKleiner et al. I(|2010D 
for the one-top codes. The BELGI-Cs-2Tops code is writ- 
ten in a quasi principal-axis-method PAM axis system 
which can be obtained by a rotation about the c-axis 
from the principal a xis system (PAM) by the angle 0, as 
shown in Eq. (9) of IQhashi et al. I (|2004f ). The value of 
9 in met hyl acetate is determ ined to be -0.0316385(95) 
radians (ITudorie et al. 1 1201 lh . U sing the values deter- 
mined bv iSheridan et al. I (|1980h in the PAM axis sys- 
tem, we obtain for the components of the electric dipolc 
moments in the quasi-PAM axis system, /i a = -0.008 De- 
bye and /if,= 1.641 Debye. 

4. RESULTS AND DISCUSSION 

The predicted frequencies and intensities, together 
with the energy of th e levels, have bee n implemented 
in the MADEX code (|Cernicharo 1 1201 2D . Each one of 
the five states, AA, AE (E 2 ), EA (Ei), EE (E 3 and 
E4) has been considered as an independent molecular 
species for the calculation of line intensities. For the ki- 
netic temperature of Orion levels above J = 30 could 
be significantly populated and some lines could also con- 
tribute to the Orion spectrum. The predictions from 
our spectroscopic code have, however, large uncertain- 
ties for larger J's (see above). The synthetic spectrum of 
methyl acetate was calculated assuming LTE conditions 
(due to the lack of collisional rates for this molecule) 
for a kinetic temperature of 150±20 K, a LSR veloc- 
ity of the cloud of 8 kms -1 , a line width of 3 kms -1 , 
and a column density for each state AA, EA, and AE 
of (10.4±1.0)xl0 14 cm' 2 , and of (5.2±1.0) x 10 14 cm~ 2 
for the Ei and E 2 states. Size (15") and offset (7") from 
the pointin g position (IRc2) o f the compact ridge com- 
ponent (see lFavre et al. 1 120111 iGenzel fc Stutzki 1119891 
iBlake et al. Ill987l ) are taken into account in our model. 
Beam dilution is corrected for each line depending on 
their frequency. Consequently, the total column density 
of methyl acetate in Orion is 4.2±0.5xl0 15 cm~ 2 which 
includes a correction factor to the partition function com- 
puted for T K =150 K and J<30 of 2.6. It includes the 
torsional excited states (estimated energies) and rota- 
tional levels up to J=65. Figures 1, 2, and 3 show se- 
lected lines from the best fit synthetic spectrum to the 
Orion data. The full list of detected lines is provided in 
Table 1 (electronic version) where we list 215 unblended 
lines as well as 163 lines moderately blended with other 
species. No unblended lines of CH 3 COOCH 3 are miss- 
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ing. The detection is fully secure taking into account 
that the systematic pattern of the lines arising from the 
different states is always present. Lines from v t = \ could 
be also detectable as the vibrational partition function 
indicates that the population of the v t =l levels from the 
two methyl groups will be a factor 1.5-2 below that of 
the ^t=0 states. Unfortunately, no accurate predictions 
are available for the torsionally excited states. 

Ethyl formate has two conformers, trans (also called 
anti) and gauche with the latter being 65±21 cm -1 
(94 K) above in energy (|Riveros fe Wilson Ill967[ ). The 
gauche conformer has two possible orientations of the 
terminal methyl group and, hence, it could be twice as 
abundant as t he trans conformer i f the energy differ- 
ence were zero (jBelloche et al. 1120091 ) . We have used the 
spectr o scopic laboratory data f rom iRiveros fe Wilson I 
(1967), iMever fc Wilson I (|1970l ). and iDemaison et al. I 
( 1984) to derive the rotational constants which were in- 
corporated into the MADEX code. The predictions were 
chec ked against the corre sponding entry in the JPL cat- 
alog {Pickett et al. 199_1). The t rans conformer of ethy l 
formate was detected in SgrB2 bv lBelloche et alTl (2009) 
but only upper limits were obtained for the gauche one. 
We have identified both conformers in the line survey of 
Orion trhough 90 features free of blending (52 trans and 
38 gauche). Selected identified lines are shown in Fig. 
4. We found that the gauche conformer has less intense 
lines than the trans one for transitions at 3 and 2 mm. 
However, at higher frequencies there are many multiplet 
transitions of the gauche that coincide in frequency and 
they become prominent in the Orion data. The best fit 
to the data (assuming the same physical conditions and 
source parameters than those of methyl acetate corre- 
sponding to the compact ridge component of Orion KL) 
provide a rotational temperature of 150±20 K, and a col- 
umn density for each conformer of (4.5±1.0) x 10 14 cm~ 2 , 
i.e., a total column density for ethyl formate ~5 times 
below that of methyl acetate. From the observed abun- 
dance ratio between the trans and gauche conformers, 
N gauche /N trans -1= 2xe~ 94/T , we derive a kinetic tem- 
perature for the emitting gas of 135±30 K. If the low- 
est energy conformer of ethyl formate, the trans one, 
was formed on the ices before the warm phase of the 
cloud, isomerization to the gauche form will require a 
high kinetic energy and, perhaps, a lo ng time to overpass 
the b arrier to isomerization of 550 K (|Riveros fe Wilson! 
119671 ). Our result points towards a fast equilibrium be- 
tween both conformers at 150 K and in a time compa- 
rable to the duration of the present warm phase of the 
cloud. If the molecule is formed in gas phase the en- 
ergy liberated in the process could help in the isomer- 
ization process and the observed conformer temperature 
will reflect the kinetic temperature of the gas. However, 
no chemical paths are included in the present chemical 
models to form neither ethyl formate nor methyl acetate. 
New laboratory experiments are needed to understand 
the way these species could be formed in gas phase or in 
ices and how conformers of the same species can equili- 
brate at a temperature close to the kinetic temperature 
of the gas. 

Methyl acetate is probably formed via multiple reac- 
tion pathways from species detected in hot cores such as 
methyl formate, acetic acid, and methanol. In icy grain 



mantles, methanol is one of the common molecular com- 
ponents. Acetic acid is a relatively large molecule that 
would be difficult to detect in the ice by IR observa- 
tions and that could be the precursor of methyl acetate 
and ethyl formate. Only formic acid (HCOOH) has been 
proposed as a possible carrier of the 7.24 /jm band to- 
ward high-mass protostars, while the 7.41 /xm band could 
be due to the formate io n (HCOO~) and acet aldehyde 
(CH 3 CHO), according to lSchutte et al. I (|1999f ). Forma- 
tion of HCOOH on a surface occurs experimentally at low 
temperatures, mainly through hydr ogenation of the HO- 
CO complex (jloppolo et al. Il201ll ). Other possi ble for- 
matio n routes in the ice are via precursor cations (|Woon I 
120111 ) or by reactions of super thermal Q( 3 P) atoms an d 
CH4 with an over coat of CO ([Madzunk ov et al. Il2010f) . 
Also photon or electron irradiation of H2 0:CO ice mix- 
tures leads to formation of formic acid among other 
products , including methanol in the case of photopro - 
cessing (jWatanabe et al. I l2007t iBennett et all 120 111 ). 
HCOOH was found to spontaneously deprotonate when 
sufficient water is present to stabilize charge transfer 
complexes. Both ammonia and water can serve as pro- 
ton acceptors, yielding am monium (NHj") and hydro- 
nium (H 3 0+) counterions (|Park fc Woon I I2006D . The 
so-formed formate ion (HCOO~) might intervene in the 
formation of species like methyl acetate in the ice ma- 
trix, but t his was to our knowled ge not confirmed exper- 
imentally. iBrouillet et al. I (|2013l ) have recently observed 
with high angular resolution CH3OCH3 and CH3OCOH 
towads Orion-IRc2 and conclude that the similarity in 
the spatial distribution of both species points towards a 
common precursor. The observation of a similar abun- 
dance for the two conformers of ethyl formate points 
to a gas phase production path rather than to a low 
temperature ice for mation mechanism . Radicals such 
as methoxy rCH.O. ICernicharo 1120121 ) or CH 3 CO could 
play an important role in the gas phase chemistry. How- 
ever, m ethoxy has been observe d only towards cold dark 
clouds (jCernicharo et al71l2012D and it is not detected in 
our line survey of Orion; CH3CO has not been detected 
yet in the ISM. 

The two conformers of ethyl formate and methyl ac- 
etate are isomers of the C3H6O2. Three additional iso- 
mers, propanoic (propionic) acid (CH3CH2COOH), hy- 
droxyacetone (CH3COCH2OH), and methoxyacetalde- 
hyde (CH3OCH2OH), could be also present in Orion. 
The three species are implemented in MADEX. For 
hydroxyacetone the available spectroscopic data have 
been summarized by iBraakman et al. 1 (|2010f) and cover 
frequencies up to 431.8 GHz (dipole moments from 
iKattiia-Ari fc Harmony I I1980D . We obtain an upper 
limit to its column density of 8xl0 13 cm~ 2 (see also 
lApponi et aT~ll2006l for an upper limit to its column den- 
sity towards SgrB2). For propionic acid and methoxyac- 
ctaldehydc we also obtain upper limits to their column 
densities of 1.6 xlO 14 and 2xl0 14 cm -2 , respectively. 
We note, however, that frequency predictions above 40 
GHz for these two molecular species are rather uncertain 
( Stifv ater WWf^t, lOuvang fc Howard Il2008t iHirano et all 
1987). Hence, of the known possible non cyclic isomers 
of C3Hg02, methyl acetate appears to be the most abun- 
dant one. Laboratory spectroscopic data are needed for 
propionic acid and methoxyacetaldehyde in order to draw 
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further conclusions on their contribution to the ice man- 
tle and gas phase chemistry of hot cores, and to the forest 
of still unknown spectral features in Orion. 
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Table 1 

Detected lines of CH3COOCH3 



J 


K a 


K c p 


J' 


K 




State 


Predicted 
freq. (MHz) 


Error 
(MHz) 


E u 
(K) 


Sij fi 2 


Observed 
freq. (MHz) 


(K) 


Blend 


13 

13 
13 
13 







1 


13 

13 
13 
13 


12 

12 
12 
12 


1 

1 
1 




12 

12 
12 
12 


E3 

E4 
EA 
E4 


82779.635 

82779.725 
82780.026 
82808.862 


.012 

.012 
.012 
.012 


28.6 

28.6 
28.6 
28.6 


31.30 

31.30 
31.30 
31.30 


82780.2 

t 
t 

82809.2 


0.03 
0.04 


CH 3 CH 2 CN 

33 so 2 

HCOOCH3 


13 
13 
13 
13 


1 
1 





13 
13 
13 

13 1 


12 
12 
12 
12 




1 
1 


12 
12 
12 

12 1 


E3 
EA 
AE 
AA 


82808.881 
82809.267 
82823.351 
82823.702 


.012 
.012 
.012 
.012 


28.6 
28.6 
28.6 
28.6 


31.30 
31.30 
31.30 
31.30 


t 
t 

82823.2 
t 


0.04 


"t = l 



Continue in the electronic version. 

Note.- Emission lines of CH3COOCH3 present in the spectral scan of Orion KL from the IRAM 30-m radio-telescope. Column 1-8 indicate 
the line transition, Col. 9 the state of the molecule, Col. 10 gives the predicted frequency in the laboratory, Col. 11 uncertainty of frequency 
predictions, Col. 12 upper level energy, Col. 13 line strength, Col. 14 observed frequency assuming a vlsr of 8 km s _1 , Col. 15 mean beam 
temperature, and Col. 16 blends, 
f Blended with previous line. 



